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Abstract. Using results on inverse spectral problems, in particular 
the so-called new wave invariants attached to a classical equilibrium, we 
show that it is possible to determine the Morse index of height func- 
tions. For compact Riemannian surfaces M C R 3 this imply that we 
can retrieve the topology (via the genus). 

Our results are independent from the choice of a metric on M and can 
be obtained from the choice of a 'generic' height-function. For surfaces 
of genus zero, diffeomorphic to a 2-sphere, the method allows to detect 
the convexity, or the local convexity of the surface. 

keywords : Micro-local analysis; Schrodinger operators; Inverse spectral 
problems. 

1. Introduction. 

1.1. Basic definitions and setting. We are here interested in applications 
of the inverse spectral problem for certain special differential operators in 
the semi-classical regime. Let X be compact (boundaryless) Riemannian 
manifold equipped with a strictly positive density and Ax the (positive) 
induced Laplace-Beltrami operator. In particular, we consider the so called 
/i-quantized Schrodinger operator: 

P h = h 2 A x + V, onL 2 (X), 

also called semi-classical Schrodinger operator, where the potential V is 
measurable and bounded from below on X. 

Remark 1. With some mild conditions on V, we could also assume that X 
is non-compact. But to simplify we stay in the compact situation. 

By a standard result, see [3], when V is bounded from below Ph has a self- 
adjoint realization on a dense subset of L 2 (X). To this quantum operator 
Ph we can associate a classical counterpart with the Hamiltonian, or total 
energy , function on the phase space: 

P(x,0 = \ml + V(x) on T*X. 

Here the notation: 

ij 
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designs the norm (or scalar product at x) induced by the Riemannian metric 
G = g^(x) at x. We note $t the Hamiltonian flow of H p = d^p.d x — d x p.d^. 

Here we are mainly interested in an asymptotic relation between the semi- 
classical eigenvalues {Xj(h) : j £ N} of P}{. 

Ph(j)j(x,h) = Xj(h)4>j(x,h), 4>j G L 2 (X), as h — > + , 

and the set of fixed point V (see below) for the map <&t (viewed as a map 
on Mt x T*X). We refer to the introduction of [7J for a general presentation 
of this kind of relation between quantum and classical mechanics. 

In the last section of this article we will work with compact orientable 
surfaces M of M 3 and Schrodinger operators Ii 2 Am + z attached to height- 
functions z on M. 

Spectral statistics. 

Consider an interval I = £2] with E\ < E2 and /(e) = [E\ — e, E2 + e\. 
For each e > the pullback p~ 1 (I(e)) is a compact subset of T*X and by 
a standard argument, see [3j, it follows that the spectrum a(Ph) H /(e) is 
discrete and consists for each h in a sequence: 

Xx(h) < X 2 (h) < ... < Xjih), 

of eigenvalues of finite multiplicities, if e and h are positive. In general 
no formula is known to compute the eigenvalues Xj (h) and to get pertinent 
information about the spectrum (and the classical dynamics) it is interesting 
to study the following spectral distributions: 

(1) Y(E,h,(p) = £ 

where ip £ 5(M) is a test function, conveniently chosen, see below. 

This kind of problem leads to a mathematically rigorous version of the so- 
called Gutzwiller formula (see |13|). a formula intensively used in physics and 
quantum-chemistry. For example see [14] for various applications in physics 
and quantum chaos with also many references. The so-called semi-classical 
approximation consists in studying the asymptotic behavior of Eq.([T]) as 
the parameter h tends to zero. For a rigorous mathematical study of this 
problem a non-exhaustive list of references is [U [U \T7\ [18] . 

Wave and new wave invariants. 

To study Eq. ([1]) , a classical approach (see [H [5] and section 3 of [7] ) is to 
study the asymptotic behavior, as h — > + , of the localized trace: 



u(E,h,t) = Tr (e(P h )e-^ Ph - E A ,G G C, 



-<oc 

y • 



I will follow now the terminology used in [15]. For E regular it is known 
since a while that £1 admits an asymptotic expansion of the form: 

00 

uj{E,h,t)~ aj(E,t)h j , as h -> + . 
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Definition 2. The coefficients a,j(E,t) are some distributions on the line 
Rt and are called wave invariants of Ph. 

When E = E c is critical certain new coefficients appear in the asymptotic 
expansion of u(E,h,t) as h tends to 0. In particular, one can predict a 
general expansion in the form: 

n— 1 oo 

(2) ui(E c ,h,t) Hk(E c ,t)hi\og(h) k , as h -> 0+ 

k=0 j=—no 

for some p £ N* , see [1] for the singularity near t = and [7J for many 
examples of new wave invariants. In this work we only use the top-order 
coefficients of Eq. @ near a non-degenerate singularity (see Theorem [S] 
below). 

Definition 3. These extra distributional coefficients appearing in Eq.^j are 
called new wave invariants. 

In general, the top order coefficient w.r.t. h of the expansion involving the 
new wave invariants contains many information on the shape of the symbol. 

2. Hypotheses and semi-classical results. 

We recall first the general result determining the wave invariants at a 
regular energy level. Consider X a closed smooth Riemannian manifold, 
n = dim(X) and V G C°°(X) a positive potential. Let {\j(h) : j G N} be 
the spectrum of the Schrodinger operator Sh = h 2 Ax + V where Ax is the 
Laplacian on X (here given as a positive operator). This spectrum of Sh is 
always discrete, and with finite multiplicities, when X is compact or if V 
is 'confining', i.e. V(x) —> oo when d(x,Xo) — > oo for some xq £ X. The 
justification is that both conditions insure that the level-sets T,e defined 
below are compact and that the resolvent of Sh is a compact operator. 

Let $ t be the Hamilton flow of H(x,£) = \ |£| \ 2 X + V(x) on T*X and given 
E > define the energy surface: 

E E = {(x,OeT*X:H(x,S) = E}. 

Recall that the flow can be viewed as a map &t '■ — > (conservation 
of the energy) . Also when the surfaces £ e are compact the general theory 
of differential equations insures that the flow is complete (property of the 
maximal solutions of a Cauchy-problem). To simplify notations we write 
z = (x, £) G T*X. We recall that E is regular (or non-critical) when dH ^ 
everywhere on Eg and critical otherwise. Below, we use the subscript E c to 
distinguish out critical values of H. The so-called period manifold of on 
Y.E is: 

V = {(T, z)£lxS £ : $ T (z) = z}. 
At a non-critical energy level. 

When E is non-critical, we have, see [TT] or [5], the following general result 
concerning the wave invariants at the energy E: 
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Theorem 4 (Semi-classical trace formula at a regular level.). 

Assume that E is regular and that the restriction of &t to Tie is a clean flow 
(see section 2 of [H]j. Then there exists a sequence of distributions on the 
real line, {^k}k, such that for every test function ip with Fourier-transform 

j=0 j=l 

Moreover, the supports of the distributions jj are contained in the sets of 
periods of the closed trajectories of <3?i onY^E- 

For a better description of the coefficients appearing in Theorem H] we 
refer to pTJ [5] (we do not need their explicit expressions here), see also 
[9J for a similar high-energy result concerning elliptic operators on a com- 
pact manifold. Also, under our hypotheses the trace in Theorem H] and the 
functional T(E,h,cp) are equal modulo a function of fast decay w.r.t. h. 
Such a coefficient, of order 0(/i°°), is negligible in semi-classical asymptotic 
expansions. 

The idea we want to use here is that by a clever choice of supp((^) we can 
eliminate the wave invariants appearing in Theorem UJ 

• If ip is flat at the origin the set: 

{{0} x T E } C V, 

does not contribute. 

• If supp((^) C [—To, To], for To small enough then no periodic orbit: 

{(T,z) :<S> T (z)=z}cV, 
will contribute to the asymptotic expansion. 
At a critical energy level. We allow now the presence of critical points 
for H and we impose the type of singularity : 

(^4i) The potential V is a Morse function on X. 

A fortiori, in / there is finitely many critical values E],, E l c and in p _1 (7) 
finitely many fixed points Zq,...,z of the energy function p. 

Remark 5. The number of critical points Zq = (xq, 0) is equal to the number 
of critical energy levels. Otherwise V would not be a Morse function on X. 

Next, we impose two conditions on our test function ip: 

(A 2 ) (p is flat at 0, i.e. c^(0) = 0, Vj 6 N. 

(.A3) For some sufficiently small T we have supp((/3) C [— T, T]. 

A fundamental property is that the singularity of T(s,h,ip) as s — > E c 
describes partially the singularity of V. In fact with conditions (A2) and 
(A3) we will only see the new wave invariants attached to the critical point 
z 3 in TtEj- We have: 
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Theorem 6 (New wave invariants at a critical level). 

Under the conditions (Ai), (A2) and (A3) we have: 

' 3=0 
The leading coefficient is of the form: 

(2tt) J |det(d$ t (^)-Id)|3 

We refer to [TTJ [TJ El [16] for a proof. Observe that (.A3) implicitly insures 
that det(d&t(zl) — Id) 7^ 0. Because of our implicit choice for tp, not all 
the new wave invariants are present in this formula. The other new wave 
invariants are studied: 

• In [4], near t = 0. 

• In [6j [16] , near a period of d$t(^). 

The explicit determination of all wave invariants, near a critical point of 
arbitrary signature, is a somehow complicated analytic problem involving 
oscillatory integrals with degenerate phases. For an operator which is not 
a Schrodinger operator some new terms can generally appear at a period of 
d$ t (z Q ) (see [6]). 

New wave invariants. In our setting, the top-order coefficient, given by 
the Duistermaat-Guillemin-Uribe density, is indeed a smooth function as 
long as we stay away from any period of the linearized flow at the point zl. 
When X = M n an explicit computation, done in [16] in a suitable system of 
linear coordinates, shows that: 

(3) dv t (zo) = — l — . 

I [] sinh(a j(z )t) H sin(ay(z )i)| 

j=l j= r +l 

We must simply retain that the density dvt(zo) determines: 

• The signature (n — r, r) of the Hessian of V at zq. 

• Eigenvalues ctj(zo). 

The last affirmation follows via Taylor-series and evaluation at several times. 

Remark 7. In general, if the metric and the height function are unknown, 
the spectral expectation determines only the numbers atj and not the respec- 
tive eigenvalues of G(xq) and d 2 V(xo). A similar indetermination is already 
valid for linear combinations of harmonic oscillators on W l . 

It follows that, when the potential is Morse- function on X, we can retrieve 
the morse index of X by several successive applications of Theorem 6: we 
have only to cross finitely many critical energy levels and to collect the index 
at each energy. 
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3. Application to surfaces. 

There is a nice application to compact smooth surfaces McM 3 equipped 
with a Riemannian metric (not necessarily the metric of M 3 restricted to 
M). We assume that M C M 3 is smooth, boundaryless, orientable and 
that M carries a smooth Riemannian metric G, fixed once for all. We take 
A m as the Laplace-Beltrami operator attached to this metric (following the 
convention of geometers we may assume that Am is positive). Let us chose 
as potential V a height function. We can assume V to be positive, this can 
always be achieved via a translation, M being compact. If we embed M in 
M 3 , via some coordinates (x, y, z), we can chose V as the projection on the z 
axis. It is a standard result of topology, see chapter 6 of [2j, that for almost 
embedding V will be a Morse function. 

Then, for any choice of a smooth Riemannian metric on M, we have: 

Proposition 8. Under the previous conditions on V, the semi-classical spec- 
trum of Ph = h 2 Apf + V(x), defined as an unbounded operator on L 2 (M), 
determines the topology of M . 

Remark 9. Observe that the knowledge of the metric is not required. We 
only need a kinetic energy operator which is micro-locally elliptic and with 
a principal symbol nowhere degenerated (see below). The knowledge of V 
is also not required. We only need to recover the number of critical points 
of V and their signature to conclude. 

Proof of Proposition [8j We will use a variational argument w.r.t. the 
energy E. Since M is compact our potential has a maximum E m3iX and it 
will be sufficient to perform spectral estimates below E m3iX . Let \j(h) be the 
spectrum of Ph, each eigenvalue being repeated according to it's multiplicity. 

Since our potential is a Morse function, by Sard's theorem, we obtain 
that the energy function p(x,£) has only finitely many critical values E 3 C , 
j £ {1, ...,N}, attached to single critical points. When supp(c/3) is small 
enough and does not contains the origin we have: 



j(E,h,tp) 



0(h°°), for E non-critical, 
cj(<£) + O(h), for E = E{ critical. 



Hence, the semi-classical spectrum determines each critical value E 3 C of V. 

Now for j fixed we can use a simple micro-local argumentation. The only 
critical point on T, E j is of the form = (x J , 0) with V(xq) = E c - We pick a 
function in tp £ Cg°(T*M) such that < ip < 1 everywhere and ip = 1 in a 
neighborhood of z J . Always with our conditions (^2) and (^3) on supp(c/3), 
we have: 

Y(E, h, if) = Tr (^ w (x,hD x M Ph ~ Ei )) +0(h°°). 
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The important fact here is that on supp(l — ip) there is no critical point of 
p. Hence with condition (A2) and (.A3) we have: 

Tr ((l-^( x ,hD x ))cp(^^)j =0(h°°), 

which easily follows from a non-stationary phase argument. On supp(^), 
which can be chosen arbitrary small up to an error of order 0(h°°), we can 
use local coordinates around x J and the Laplace operator has the form: 

h 2S P r ® 1 ■■ ® +V 
^ri dxi y/g^ dxj 

= - h2 E^irir + v + h2 Y; ^7= 9 gi ^ x) ^ 

Here we write the metric G = g^ , G^ 1 = gij and g = detG. Hence, in the 
sense of the /i-calculus, we have ph := po + hp\ with a principal symbol: 

pofaQ = J2gij( x )£j£i + v(x), 

and a sub-principal symbol: 

Observe that p\ = at every point where £ = and the sub-principal symbol 
pi will play no essential role for the estimates below. You can also use the 
convention that h 2 A is the quantization of ||£|| 2 . This makes no difference 
for the spectral estimates below. 

Next, since V is independent of the choice of the metric on M we can 
freely assume Q that near the origin: 

9ij (x) = ld + 0(\\x\\). 

In this system of local coordinates we obtain that: 

Tr L^( x ,hD x M^^) \ ~Cg(0) + 0(/»), 



where: 



C 3 o (0) = C J ut(4)<p(t)dt, C G C* 



*Of course this choice is not isometric, in reality here should be the positive eigenvalues 
of G at Xq. But this choice is clearly sufficient to attain our objective. Observe that one 
could also locally trivialize the metric. This method is a well known trick to derive Morse 
inequalities, e.g. for Witten-Laplacians. 
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In our setting we can apply standard results on the linearized flow, see pQ, 
to compute the density vt(z J ). This density, of the form given by Eq.([3|), 
determines the number of positive and negative eigenvalues of V at Xq. 

In general position, define the Morse index of a critical point xq as the 
dimension of the negative eigenspace in xq. If z is Morse function we denote 
by Nj(z) the number of critical points of z with index j. Since for a surface 
we have only j = 0, 1, 2 we can retrieve the Euler characteristic: 

N (z) - N^z) + N 2 {z) = X (M) = 2(1 - g(M)), 

where g{M) is the genus of M. ■ 

Remark 10. Of course this approach is still valid in dimension n > 2 
but then the genus is no more a sufficient topological invariant. For a nice 
overview on Morse theory and indexes see [2] or [10] for surfaces. Observe 
that, for an unknown metric and an unknown height function, we can still 
retrieve the Morse index of V but not the Hessian of V . This is because in 
formula [3] we can only retrieve the ratios ctj(zo) of eigenvalues of G(xq) and 
c1 2 V(xq) in a given system of coordinates. 

The notion of Morse-Smale function (see, e.g., [2] p. 158) is here cen- 
tral. Morse-Smale functions are moreover dense in every C r -spaces (r > 1) 
(Kupka-Smale-Theorem, p. 159 and remark 6.7 p. 160 in |2j). To a Smale- 
Morse function is attached the Morse-Smale-Floer complex and this com- 
plex is isomorphic to the complex giving the singular homology (Theorem 
on Morse homology, Theorem 7.4 of [2]). For a surface it is known that the 
topology is given by the genus or the Euler characteristic and this one is 
also given by the Euler- Poincare characteristic of the complex of homology. 
That a certain Morse-function determines the topology of a surface is also 
contained in the book [10] page 70. 

About the choice of a height function. 

For certain simple surfaces, e.g., convex surfaces diffeomorphic to a 2-sphere, 
there is no bad embedding since the height-function z always shows up a 
strict minimum and a strict maximum. Such a height function is also a 
perfect Morse function, i.e. a Morse function with exactly 2 critical points 
attached respectively to a strict minimum and a strict maximum. For a sur- 
face of genus we can still have several critical points at the same critical 
value. This is not generic and unstable under a small perturbation of z. 

For surfaces of higher genus some 'bad' embeddings are possible if the 
height function is chosen transverse to a level set (non-generic choice). 

The first embedding (here z is the axe of symmetry) is not favorable: the 
set of critical point consists of 2 circles. These circles are manifolds of critical 
points attached respectively to a maxima and a minima of V of energies 
E m m < -Emax- I n that situation, a compact manifold of critical points of 
dimension 1, we can here anyhow apply the results of [4] or [16]. If we 
still assume that conditions (^2) and (^3) are satisfied, each of these circles 
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Figure 1. Embedding of a torus in R 3 . 



contributes as a 1-dimensional submanifold in the 4-dimensional phase space: 



for the circle of minima. Observe that the order w.r.t. h is now —1/2. 

Remark 11. Both formulae for Y(£ max , .) and T(E m [ n , .) easily follow from 
an application of a stationary phase method with a compact manifold of 
critical point and a non-degenerate transverse Hessian. 

Here a{9) is the ratio of the eigenvalues of the linearized operator in the 
transverse direction to Er . ~ ~ S 1 evaluated at the point € § . 

Observe that a(0) is negative for E = i? mffl (hyperbolic flow) and positive 
for E = -Emin (periodic flow). Observe that, only from the spectral estimates, 
we can still see: 

• Hyperbolic contributions: unstable equilibria at the maximal energy. 

• Trigonometric contributions: stable equilibria at the minimal energy. 

The previous situation can be generalized for a smooth curve 7, necessarily 
isomorphic to S , of critical points with a non-degenerate transverse Hessian 
at each point of 7. 

Remark 12. For a generic choice of the metric G on T 2 the function 
9 1 — y a(9) is not constant along § . Unfortunately, the spectrum of the 
associated Schrodinger operator —h 2 A g + z only determines the average of 
the density along the circles. To get a better description here requires to 
perform eigenfunction estimates. See, e.g. j4] for this point. 




tesupp(^) 



for the circle of minima and: 




egs 1 tesupp^) 
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For the second embedding, where x is the axe of symmetry, z is a Morse 
function and we meet successively the critical points: 

• a singularity of type (0, 2) : strict minimum, 

• a singularity of type (1,1) : first saddle point, 

• a singularity of type (1,1) : second saddle point, 

• a singularity of type (2, 0) : strict maximum. 
This gives: 

X (T 2 ) = l- 2 + l = 0=> g{T 2 ) = 1. 
The situation of the second embedding is generic and stable (e.g., w.r.t. a 
little deformation of the height function) . This allows to retrieve T 2 , up to 
a smooth deformation. 

Remark 13. Results concerning Morse-functions are not specific to height- 
functions and Proposition [8] can be generalized to /i 2 Ajf + V where V is a 
Morse-function on M. The interest here is the evident physical interpreta- 
tion: the choice of V(x) = z(x), x G M is equivalent to put a particle, forced 
to move freely on M along geodesies, in a constant gravitation field. Also a 
height-function gives a function independent of the choice of the metric on 
M. For example, this is not the case of a potential V(x) depending (locally) 
on the geodesic distance d(xo, x) on M. Such a potential depends on G and 
can be singular at conjugate points. 

Convexity and measure. 

Assume that g(M) = then if M is convex every choice of a height-function 
gives a perfect Morse function. But if M is not convex certain choice of the 
potential give locally a number of critical points greater than 2, with the 
same Morse index. 




Figure 2. A non-convex surface: locally we see 3 critical points. 

Observe that we still get 3 critical points by a small change for z and that 
we get an open condition. In general to detect the non-convexity of M does 
not require global spectral-estimates. Finally, for any gender g{M) observe 
that global spectral-estimates provide a lower bound for the euclidian di- 
ameter of M since we can retrieve the maximum and the minimum of the 
restriction z\m- This is just a lower bound as shows figure 1. 
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From the point of view of statistical mechanics it could be interesting to 
put a probability measure [x(z) on all choice possible for z and to average 
the spectral estimates with respect to /x. This is here simply equivalent to 
chose a probability measure on S 2 since the full problem is invariant under 
translation. A similar construction is possible for the choice of the metric 
G: if G is in a bounded set of metrics G a , estimates given by conditions 
(A2) and (.A3) are still globally valid and so are our conclusions. One could 
average the results with respect to some probability measure fJ,(a). 

We could also obtain the contributions of a surface M carrying a flat section 
in the following sense: 

There exist an open subset [/CM 3 and a two-dimensional 
plane PcK 3 such that U DM = U DP. 

The flat section can be interpreted a 2-dimensional subset of critical points 
when the height function z is chosen transversally to this section. The 
associated result is simply the Lebesgue measure of the flat section. Observe 
that, as predicted by the general theory of Morse- functions, this situation is 
not generic and not stable under a small perturbation of z. 

Final Remark. At a first look it might seem childish to use semi-classical 
methods. But the 'high-energy' method (see e.g. [9]) is not working: when 
the energy E is larger than the maximum of the potential E max we have 
that the kinetic energy is bounded from below by \\^\\ x > E — E max > 0. 
By ellipticity of the Laplacian, it is not possible to produce any new wave 
invariant in this regime. 
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